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مدل‌سازی پارامترهای جریان در مسیر تنداب سرریزهای با تغیبر شیب ناگهانی 
(مطالعة موردی: سرریز سد سور کت) * 
«یادداشت پژوهشی؛ 
سهیل عباسی" حسین صمدی بروجنی" ‏ روح‌الله فتاحی نافچی"" بهزاد قربانی"" 

چکیده یکی از موارد مهم در طراحی تنداب سرریزهاء بررسی حصوصیات و شرایط جریان در مسیر تنداب به‌واسطهُ فوق بحرانی بودن 
جریان است. در تحفیق حاضر ویژگی‌های هیدرولیکی جریان در تنداب سرریز سد سورک بااستفاده از مدل‌سازی فیزیکی و عددی مورد 
بررسی قرار گرفته‌است. مدل از جنس پلکسی‌گلاس با مقیاس ۱:۵۰ ساحته شد و برای مدل‌سازی ریاضی از نرم‌افزار فلوئنت استفاده گردید. 
آزمایش‌ها در ۵ دبی متفاوت انجام گرفت و پارابترهای فشار, سرعت و عمق جریان در محور مرکزی و کنار دیواره برداشت گردید. نتایج 
نشان داد در بین پارامترهای اندازه‌گیری‌شده نرم‌افزار فلوگنت» پارامتر سرعت با کمترین درصد نعطا ( 8/۳ درصد) و پارامتر عمق با بيشترین 
درصد عطا (۱۲/۹ درصد) شبیه‌سازی شده‌است. هم‌چنین ضریب نش- سائکلیف و ضرپب تعیین محاسبه‌شده برای پارامترهای مورد نظر 
نشان دادند که نرم‌فزار فلوئنت توانایی شبیه‌سازی جریان را برروی سرریزهای با تغییر شیب ناگهانی دارا می‌باشد. بانوجه به نتایج شبیه- 
سازی‌شده و معیارهای 10/51 و ۱۷5 می‌نوان دریافت در نرم‌افزار فلوئنت درصورت استفاده از روش حل عددی حجم محدود؛ مدل 
آشفتگی ٩-۵‏ 1۷0 روش تصحیح فشار 3114۳1 و روش دو فازی ۷0۳ بیشترین دقت در محاسبه پارامترهای فشار» سرعت و عمق جریان 
آب در تنداب سرریز مورد مطالعه حاصل می‌گردد. 


واژه‌های کلیدی فشا سرعت. سرریز. مدل فیزیکی. مدل‌سازی عددی. 
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۱3۹۹ مدلسازی پارامترهای جریان در مسیر تنداب سرریز های با تغییر شیب اگهانی.. 


مقد مه 

سرریز یکی از سازه‌های مهم هر سد برای عبور 
سیلاب‌ها می‌باشد که تنداب آن باتوجه به شرایط 
توپوگرافی و شیب طبیعی زمین در مسیر ممکن است با 
شیب‌های مختلف طراحی و اجرا شود این امر تغییر 
شیب‌هایی ناگهانی در تنداب به‌وجود می‌آورد. بدیهی 
است در این وضعیت. سازهٌ تنداب باید از نظر شرایط 
هیدرولیکی. به‌گونه‌ای طراحی شود که پدیده‌هایی 
همچون کاویتاسیون. ایجاد امواج ثانویه و.. در محل 
تغییر شیب ناگهانی حادث نشوند و شرایط عبور ایمن 
سیلاب فراهم باشد [1]. 

براساس اطلاعات موجود. در موارد زیادی دلیل 
تخخریب سرریزهاء طراحی نامناسب تنداب‌ها بوده‌است. 
به دلیل اقتصادی و هیدرولیکی, مناسب‌تر است که 
تنداب. مستقیم و شیب کف آن از شیب بحرانی بیشتر 
باشد. با این حال شرایط زمین‌شناسی و مورفولوژیکی 
شیب دره و تکیه‌گاه که تنداب باید از آن پیروی کند 
گاهی ممکن است موجب شود که در مسیر تنداب حم 
حاصل شود و هم‌چنین شیب کف آن در مسیر تغیسر 
نماید. به‌طور کلی اصول زير در طراحی تنداب مورد 
توجه قرار می‌گیرند [2]. 
. هرگونه هدایت و یا تغییر مسیر جریان باید محدود 
به مکان‌هایی واقع در بالادست سازه کنترل (جایی که 
سرعت جریان کم است) باشد. 
۲ در منطقه‌ای که جریان شتاب می‌گیرد. هدایت جریان 
باید مستقیم و متقارن باشد. 
۳ هنگام عبور جریان از روی ساز؛ کنترل جریان 
بهسرعت شتاب می‌گیرد. پایین‌دست آستانه. جریان فوق 
بحرانی است» سرعت آن افزایش می‌یابد و امکان ایجاد 
کاویتاسیون را فراهم می‌کند؛ بنابراین هرگونه تغییر در 
انحنای دیوار هادی جریان و کف تنداب واقع در 
پایین‌دست آستانه» باید با دقت زیاد انجام شود. 
ازآن‌جایی که پهنه‌های بالادست آستانه وابسته به دبی 
عبوری از مقطع کنترل می‌باشد و اغلب پهنا در مقطع 
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کنترل نسبت‌به تنداب بزرگ‌تر است. بنابراین یک تبدیل 
بین این دو لازم است؛ البته پهنای تنداب به‌وسيپلة 
عواملی از قبیل توپوگرافی هزينة ساخت و پهنای 
رودخانة پایاب محدود می‌شود. 

بنابراین بررسی رفتار هیدرولیکی و شناخحت 
پارامترهای هیدرولیکی سرریز در دبی‌های مختلف 
به‌کمک مدل‌های فیزیکی و عددی امکان‌پذیر است. 
تغییرشیب‌های ناگهانی در این نوع سازه‌ها می‌تواند 
زنمهشاو بق نک ضربان برد و اگوی سرنان ونر 
دهد که این امر می‌تواند مسائلی را در پایداری تنداب 
ایجاد نماید؛ بنابراین بررسی پارامترهای جریان (توزیع 
فشار و سرعت) در تغییر شیب‌های ناگهانی بهحصوص 
در جریان‌های با سرعت زیاد که می‌تواند موجب 
عملکرد نامطلوب هیدرولیک جریان و ایجاد 
خسارت‌های شدید و پا شکست سازه گردد. کاملا 
الزامی است [1]. مهری [۳] به شبیه‌سازی رفتار 
هیدرولیکی جریان برروی سرریز سد بالارود بااستفاده 
از مدل فیزیکی با مقیاس کوچک پرداخت. در این 
تحقیق مدل فیزیکی با مقیاس ۶۰ ساخته شد و در 
۶ دبی مختلف پارامترهای مهم جریان شامل عمق. 
سرعت. فشار, نوسانات لحظه‌ای فشان مشخصات جت 
پرتابی و بی‌نظمی‌های جریان با هدف حصول اطمینان 
ز عملکرد مطلوب سرریز در دوران بهربرداری مورد 
ندازه گیری و مشاهده قرار گرفت. اآهم نتایج وی نشان 
داد که تراز حداکثر بازشدگی دریچه‌های سرریز و 
رتفاع دیواره‌های جام پرتابی نیاز به اصلاح دارند و 
باتوجه به این که ضریب خوردگی اندازه‌گیری‌شده بیش 
ز مقدار بحرانی آن است احتمال به‌وجود آمدن پديدة 


کاویتاسیون وجود ندارد. عارف پور [1] به بررسی 
هیدرودینامیک جریان و منحنی تراژکتوری جت 
پایین‌دست در پرتابة جامی‌شکل سرریز سد بالارود در 
مقیاس ۱:4۰ پرداخته‌است که نتایج وی حاکی از قابل 
قبول بودن دامن تغییرات فشارهای دینامیکی بوده‌است. 


دورقی [4] به بررسی رفتار هیدرودینامیک آب برروی 
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سهیل عباسی- حسین صمدی بروجنی- روحاله فناحی افچی- بهزاد قربالی ۱۳۷ 


سرریز اوجی سه‌دهانه و تعیین رابطهٌ دبی اشل و 
بهینه‌سازی دیواره‌های هدایت سرریز سد بالارود در 
مقیاس ۱:۶۰ پرداخت. در این تحقیق وی در ابتدا به 
بهینه‌سازی دیواره‌های هدایت سرریز پرداخت و 
بهینه‌سازی هم به‌صورت کمی و هم به‌صورت کیفی 
(مشاهدة الگوی جریان. امواجم عرضی. تلاطم و 
جداشدگی جریان) انجام شد. این تحقیق در 1 دبی 
مختلف پارامترهای مهم جریان شامل عمق. سرعت؛ 
کار اشتایکن توسانات بط ای فیار او بیط ها 
جریان را با هدف حصول اطمینان از عملکرد مطلوب 
سرریز در دوران بهره‌برداری. مورد اندازه‌گیری و 
مشاهده قرار داد. بررسی عمق آب روی سرریز و 
کنایت دیا ها سفن کیک زاین ۱۷ 
ثانیه بلافاصله بعد از تبدیل ابتدایی. عمق آب به بالای 
دیوارة سرریز می‌رسد که نیازمند اصلاح می‌باشد. 
حداقل ضریب خوردگی که در آستانةٌ سرریز و در دبی 
۳ لیتر بر انیه رخ داد. ۰/۸۷ بود ولی باتوجه به 
ضریب خوردگی بحرانی (0) احتمال رخ دادن پدید؛ 
کاویتاسیون در طول محور سرریز وجود ندارد. 

علاوه‌بر مدل‌های فیزیکی, درمورد استفاده از 
مدل‌های عددی و کامپیوتری دربارة جریان روی 
سرریزها کارهای زیادی صورت گرفته‌است. ساویج و 
جانسون [6] به بررسی پارامترهای جریان برروی یک 
سرریز اوجی به‌وسیلةٌ مدل فیزیکی و مدل عددی 
(30 ۳۱۵۷) با روش متوسط‌گیری 


(وعاماه -عنسعل ۸7۵۲2260 -600[06) _ پر داختند. 


رینولدز 


آنها اظهار داشتند که همواره روش‌های عددی 
درصورتی که معادلات حل جریان به‌درستی انتخاب 
کردند جایگزین خوبی برای مدل‌های فیزیکی می‌باشند. 
تال از فنکر فیتت ای ال ضل ی «طیرا2 شوت بر 
هزینه و زمان را برشمردند. کسدی [7] براساس 
جریان پتانسیل. به‌صورت دوبعدی فشار روی تاج 
سرریز را به‌وسیلةٌ مدل عددی بررسی کرد که نتایج آن 
تطابق خوبی با نتایج آزمایشگاهی داشت. در اولسن و 
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کالویج [8] بااستفاده از معادلات رینولدز 
ميانگین گیری‌شده در زمان و نیز معادلات ع- 6 استاندارد 
به‌روش حجم محدود. جریان عبوری از روی سرریز 
را به‌صورت دوبعدی و سه‌بعدی تحلیل کردند. 
یوامی و همکاران [9] با ترکیب دو روش حجم 
محدود و المان محدود. جربان سطح آزاد روی 
سرریز را به‌صورت دوبعدی با شبکه نامنظم 
مثلثی‌شکل شبیه‌سازی کردند. چنرسری و همکاران [10] 
به بررسی جریان در یک سرریز پله‌کانی و هم‌چنین 
بدون پله بااستفاده از روش‌های عددی پرداختند. این 
محققان برای اجرای مدل عددی» روش ۲ را برای 
حل جریان دوفازی (آب و هوا) و مدل ع ->1 را برای 
حل آشفتگی جریان انتخاب کردند. حیدری [11] نیز 
جربان آشفتة عبوری روی سرریز اوجی را با 
نرم‌افزار فلوئنت به‌صورت دوبعدی بررسی کرد. 

در این راستا در این پژوهش به بررسی پارامترهای 
هیدرولیکی جریان در سرریز سد سورک و تغییر شیب 
ناگهانی موجود در تنداب این سازه با بهره‌گیری از مدل 
فیزیکی و مدل عددی فلوئنت. در ء دبی مختلف 
پرداخته شده‌است. ضمن اندازه‌گیری مشخصه‌های 
جریان در طول سرریز بررسی تغییرات فشار تغییرات 
سرعت و تغییرات عمق آب بااستفاده از شبیه‌سازی 
جریان در نرم‌افزار فلوئنت مورد تجزیه و تحلیل قرار 
گرفته و کاربرد این مدل در جریان از روی این قبیل 
سازه‌ها ارزیابی شده‌است. 


روش تحقیق 
موقعیت و مشخصات سد و سرریز سد سورک 
کر سوکان مه اکن با همه ری اک هرا 
یک سرریز اوجی با تنداب می‌باشد. این سد در نزدیکی 
روستای سورک در 1۵ کیلومتری شهرکرد. مرکز استان 
چهارمحال و بختیاری قرار دارد. حداکثر ارتفاع سد از 
بستر ۳۹ مت طول تاج سد ۰ متر» حجم کل بدنه این 
سد ۱۲۵۰۰۰۰ مترمکعب: حجم مفید مخزن ۷01 ۲۵ 
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و حجم مرده آن ۸101 ۱/۱۵ است. عرض آستانة 
سرریز و تنداب برابر ۲۰ متر و شیب وجه بالادست آن 
می‌باشد. باتوجه به مطالعات هیدرولوژی منطقه دبی 
طرح این سرریز ۲۳۱ مترمکعب برثانیه درنظر گرفته 
شدءاست. لاب این وتان ترس 
نسبت‌به افق دارد که در قسمت انتها به‌طور ناگهانی با 
زاویٌ ۲۸ درجه نسبت‌به افق به حوضچ آرامش متصل 
می گردد. 


ساخت مدل فیزیکی سرریز 
آزمایش‌های این تحقیق در فلوم شیب‌پذیر به‌طول ۲۰ 
متر و ارتفاع و عرض ۰۰ سانتی‌متر مستقر در آزمایشگاه 
هیدرولیک دانشگاه شهرکرد انجام گرفته‌است به‌منظور 
انجام آزمایش‌ها. مدل فیزیکی درانتهای اين فلوم نصب 
شده‌است. سیستم تغذیه و گردش آب نیز به تبع همان 
سیستم فلوم مذکور می‌باشد که شامل ایستگاه پمپاژ با 
ظرفیت ۷ لیتر بر ثانیه است. برای اندازه‌گیری دبی 
جریان. از سرریز مثلثی واقع در انتهای سیستم استفاده 
می‌شود. باتوجه به امکانات آزمایشگاه» مقیاس مدل 
طوری تعیین شده‌است که اولاً فضای آزمایشگاهی از 
نظر هندسی کفایت ابعاد را داشته باشده ثانیاً حداکثر 
دبی مورد نیاز در آزمایش مدل در رنج دبی سیستم 
تغذیه قرار گیرد. باتوجه به مطالب گفته‌شده و بررسی- 
های صورت گرفته مقیاس مدل ۱:۵۰ انتخاب شده‌است 
که با این مقیاس امکان تشابه چه از نظر هندسی و چه 
از نظر دینامیکی فراهم می‌گردد. در اين تحقیق برای 
ساخت مدل از مصالحی با جنس پلگسی گلاس استفاده 
شده‌است. صفحات برش‌خوردة پلگسی به هم وصل 
می‌شوند و روی شاسی فلزی قرار داده می‌شوند. برای 
اندازه‌گیری فشار» در سرتاسر سرریز در مناطق با تغییر 
یاف دره شیب رای تفای با ابقر 
نصب می‌شوند. برای سهولت کار پیزومترها روی یک 
تابلو قرار داده می‌شوند. برای رسیدن به دقت بیشتر در 
اندازه‌گیری ارتفاع آب در پیزومترهاء تابلو با زاوية ۳۰ 
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درجه نسبت‌به افق نصب می‌گردد. درنهایت با برقراری 
۶ دبی مختلف. سرعت جریان به‌وسیله لوله پیتوت. 
فشار به‌وسیلةٌ پیزومتر و پارامتر های مهم دیگر مثل 
عمق با وسایل موجود (عمق‌سنج) اندازه‌گیری می‌گردد. 
شکل (۱) طرح شماتیک مدل سرریز سد سورک را 
نشان می‌دهد. در شکل (۲) موقعیت اندازه‌گیری پارامتر- 


ها نشان داده شده‌است. 


شکل ۱ تصاویری از مدل فیزیکی سرریز سد سورک و تابلوی 


پیزومترها 


شکل ۲ مکان اندازه‌گیری پارامترها 


تحلیل ابمادی 
جریان در سرریزهای همراه تنداب براساس خصوصیات 
هندسی, سینماتیکی و دینامیکی جریان به‌صورت تابع 


(۱) به‌دست آمد: 


سال سی و یکم» شمارُ چهار ۱۳۹۷ 


سهیل عباسی- حسین صمدی بروجنی- روحالّه فتاحی افچی- بهزاد قربای ۱۳۹۹ 


)۱( 20 (2 ولا و۵ و0 ,0۵0 ٩9,‏ بط ,ظ ,ظ ,۷ ,۷) ۲ 
که در آن ظ عرض کانال. ۷سرعت جریان. ۲ 
فشار وارد روی سطح سرریزء 5 شیب طولی کانال. 0 
دبی جریان» ه کشش سطحی, ۸ لزجت دینامیکی و م 
چگالی مخصوص سیال می‌باشند. ]1 ارتفاع استاتیکی 
روی سرریز و ۷ عمق آب است. تحلیل ابعادی 
بااستفاده از تثوری 7 باکینگهام صورت می‌گیرد و 
پارامترهای موثر جریان به‌صورت معادلة (۲) به‌دست 
ی ان 
1 ظ 6 ۷ 0 


9 ,9 ببس رس 
۷ ۷۲ 0۷۰۲۷۲ ۶۷۲ للع +۷ 


9 


باتوجه به این که عرض در کانال تغییر نمی کند و 
شیب‌های موجود در مدل در طول آزمایش‌ها ثابت 
مسائل مهم در مدل‌سازی سازه‌های آبی به حداقل 
رساندن اثر مقیاس می‌باشد و همواره باید توجه نمود 
برای ناچیز کردن اثر ویسکوزیته باید عدد رینولدز 
حداقل *3۰ و برای به حداقل رساندن اثر کشش سطحی 
عدد وبر برابر ۱۰۰ باشد تا بتوان اثرات نامطلوب مقیاس 
را به حداقل رساند. گونزالس و جانسون [12] و باس و 
هاگر [13] باتوجه به نتایج این پژوهش و وجود شرایط 
ذکرشده در این پژورهش کشش سطحی و نیروی 
ویسکوزیته حداقل تأثیر را در اندازه‌گیری‌ها دارند. 


مدل عددی فلوئنت 

در تحقیق حاضر جریان برروی سرریز سد سورک 
بااستفاده از مدل فیزیکی بررسی و سپس با مدل 
فلوئنت به‌روش حجم 
شبیه‌سازی جریان در سرریز سد سورک با حل معادلات 
پیوستگی و مومنتوم بااستفاده از مدل‌های آشفتگی, 
روش‌های بررسی جریان چند فازی» روش حل عددی 
حجم محدود و الگوریتم‌های حل همزمان سرعت 
فشار 513118 انجام شد و نتایج با مدل فیزیکی مقایسه 
گردید. 


محدود شبیه‌سازی شد. 


سال سی و یکم» شمار چهار, ۱۳۹۷ 


برای شبیه‌سازی جریان آشفتهُ چند فازی از مدل 
۳ و برای شبیه‌سازی جریان آشفته از روش آشفتگی 
عع و مدل 130 استفاده می‌شود. در اين مدل حل 
معادلات مربوط به‌روش حجم محدود می‌باشد. روش 
۳ که غالباً به‌عنوان روش حجم سیال نامیده می- 
شیم بر این اسان زار است که فو باسته سال: نا 
هم ترکیب نمی‌شوند. برای هر فاز یک متغیر در مدل 
درنظر گرفته می‌شود که همان نسبت حجمی آن فاز در 
سلول محاسباتی است. در هر حجم کنترل مجموع 
تشت‌های: سسی. امین فازها سک سس ناش 
چن و همکاران [14] بیان می‌کنند که مرز بین دو فاز 
در روش حجم سیال نقطه‌ای درنظر گرفته می‌شود 
که مقدار عددی پارامتر کسر حجمی در آن ۰/۵ 
باشد. اگر نسبت حجمی سیال در سلول به‌صورت 
0 نشان داده شود یکی از سه شرط (الف) 0 2 28 
سلول از سیال خالی می‌باشد. (ب) 1 < 24 سلول از 
سیال پر می‌باشد. (ج( 0۱ 2 > ۰ سلول دارای سطح 
مشترک بین دو يا چند سیال است. وجود دارد [115. 
اندرسون و همکاران [16] از معادلات پیوستگی و 
مومنتوم ناویر- استوکس به‌عنوان معادلة حساکم بسر 
جریان در سرریز استفاده کردند. در این تحقیق از 
روش متوسط رینولدزی برای حل معادلات در 
جریان آشفته استفاده شده‌است. معادلة پیوستگی 
رینولدز به‌صورت زیر می‌باشد. 
۳ 0( + روم با 
1 ۵ ۵ 
فرم کلی معادلهُ مومنتوم ناویر- استوکس متوسط 


رینولدزی به‌صورت زیر است: 
+[ر +( جی) ٩‏ 
0 
7 ۳ و 8 
,۲ -(0 + +( +( +6۲( +ع)ات 
3 0 
سس ۲ / مه ح ۵ 
+[(۲ +6۲( 0ص 


0/۱ ۱ / ۱/۳ 15 5 
[ بت - ( + + )0 07 
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۱۷۰ فلل‌باوین با رآظرهای جریان در فسیر اب سرویرفای نا تقیی شیب نا گهانی: 


مدل 1 درواقع یکی از مدل‌های آشفتگی در حل 
معادلات پیوستگی و مومنتوم ناویر- استوکس متوسط 
رینولدزی به‌شمار می‌رود. که بر سسه 
نوع استاندارد. ۳ و ۳621122016 می‌باشد. این 
مدل دارای دو معادله براساس 1 و ع است. در این 
تلفات انرژی جنبشی اختشاش می‌باشد. چادری [17] 
معادلات انتفقال در مسدل ۱۱6 را به‌صورت زیر 


ارائه داده‌است. 


8 8 
( 1610 + (16 
(رتکام) 2 (0ام) 7 
8 
- 08 + )+ 0 متا برچ 
6 
8 8 8 
5 + (میلا,2) 3 < (,ع0) ۹ + (28) ۹ 


2 
3 
و۳ (تای )+ هب ًِ 


در معادلات مذکور ع2 و 2 به‌ترتیب عکس اعداد 
اغتشاش پرانتل برای ۶ و > و #ءلل ویسکوزيتة موثر 
هی بل 

برای معادلات مومنتوم. انرژی جنبشی اغتشاش 
و نرخ استهلاک انرژی 
گسسته‌سازی مرتبه‌دوم بالادست و برای معادلة 


اغتشاش. از روش 


کسر حجمی از روش بازسازی هندسی و برای 
معادلة فشار از روش :۳8۳8۲0 استفاهه شنده‌است: 


شبکه‌بندی مدل و شرایط مرزی 
مدل سرریز سد سورک به‌کمک نرم‌افزار گمبیت و 
بااستفاده از المان‌های مثلثی مش‌بندی شد. به‌منظور 
انجام فرایند مستقل از مش سه مدل با تعداد المان‌های 
۲ و ۷۵۶۳۲ و 1۶۱۲۶ و ۶71۳۱ اجرا شد. در 
مقايسةٌ نتایج دو مدل اول اختلاف محسوسی وجود 
نداشت. باکاهش تعداد المان‌های شبکه به 1۶7۲۶ و 
۱ سلول اختلاف نتایج مدل فیزیکی و عددی 
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بیشتر گردید. به‌طوری که در مدل با تعداد المان ۶87۳۱ 
خطا مدل به‌صورت چشم‌گیری افزایش پیدا کرد؛ 
به‌همین‌منظور برای کلیةٌ مدل‌ها شبکهٌ یکسان با تعداد 
المان ۷۵۶۳۲ درنظر گرفته شد و برای شبیه‌سازی دقیق- 
تر از المان لایه‌مرزی نیز در شبیه‌سازی عددی استفاده 
گردید. 

شرایط مرزی مدل عددی همانند شرایط مدل فیزیکی 
درنظر گرفته شد که در شکل (۳) می‌توان نمای 
شبکه‌بندی مدل و شرایط مرزی را مشاهده نمود. 


۵۲255۷۲6-1 


دروووورهب رود 2۸رمداج/رن رود 
20 ۱۱۵-2۱55 13 


شکل ۳ شبکه‌بندی و شرایط مرزی درنظر گرفته شده برای اجرای 


مدل عددی 


صحت سنجی و درصد خطای مدل فلوئنت 
برای صحت‌سنجی مدل تعدادی از نتایج آزمایش‌هایی 
که در واسنجی شرکت داده نشده بودند استفاده گردید. 
در این تحقیق برای صحت‌سنجی (تشخیص میزان 
تقارب نتایج محاسبه‌ای و مشاهده‌ای) از معیار نش- 
ساتکلیف (معادلةٌ ۷) استفاده شد؛ این شاخص در بازه 
منفی بی‌نهایت و یک تغییرپذیر است و هر چه مقدار آن 
به عدد ۱ نزدیک باشد بیانگر تقارب بیشتر است و 
مقادیر کمتر از نیم برای این شاخص نشانهُ عدم تطابق و 
تقارب دو دسته داده خواهد بود. میانگین ضریب نش 
برای داده‌های آزمایشگاهی و عددی در جدول (۱) ارائه 


نله اشتتا 


سال سی و یکم شمارة چهار, ۱۳۹۷ 


سهیل عباسی- حسین صمدی بروجنی- روح‌اله فناحی افچی- بهزاد قربالی 2 


(۷) 


در رابطة (0 8 ضریب نش» :5 : مقدار 
مشاهده‌ای» معن: مقدار محاسبه‌ای» ۲: میانگین مقدار 
مشاهده‌ای می‌باشد. 

به‌منظور محاسبة درصد خطای مدل فلوئنئت 
نسبت‌به مسدل آزمایشگاهی از رابطة (۸) استفاده 
شده‌است که نتایج آن در جدول (۱) قابل مشاهده است. 
« 


۸0 100 (۳ 


1 
در رابطة فوق ,ره مقدار پارامتر اندازه‌گیری‌شده در 
مدل عددی. 1 مقدار پارامتر محاسبه‌شده در مدل 


و ۸0 درصد خحطای مدل می‌باشد. 


جدول ۱ میانگین ضریب نش - ساتکلیف و میانگین خطا 


نام پارامتر ی میانگین درصد خطا 
ساتکلیف 
فشار ۸۵ ۳ 
سرعت ۸۷ 1/۳ 
عم ان ۹۷ ۱3/۹ 


نتایج و بحث 
اندازه گیری پارامترهای جریان در مدل فیزیکی در چهار 
دپی ۷۳/۶۱ و ۱۱۰/۸۰ و ۱۷۲/۷۳ و ۲۳۷/۷۶ مترمکعب 
بر انیه و در ۵۲ موقعیت در امتداد خحط وسط و 
کناره‌های مدل صورت گرفت مدل عددی نیز بااستفاده 
ثانیه کالیبره و بااستفاده از داده‌های مرتبط با دبی‌های 
دیگر مدل صحت سنجی شاه درادامه به بررسی نتایج 


تست املم از مدل فیزیکی و عددی پرداخته شده‌است. 


بررسی تغییرات عمق جریان 
تغییرات عمق آب از تاج سرریز تا انتهای تنداب 
به‌صورت کاهشی مشاهده 7 این روند در مدل 


سال سی و یکم» شمارة چهار, ۱۳۹۷ 


عددی نیز به همین صورت شبیه‌سازی شد. در چهار 
دبی ۲۳۶/۷۳ ۱۷۷/۷۳ ۱۱۰/۸۰ ۷۳/۶۱ مترمکعب بر 
ثانیه مقدار ارتفاع آب نسبت‌به تاج سرریز در مدل 
فیزیکی پس از تبدیل به نمونة اصلی به‌ترتیب ۲/۵۳ ۲ 
۷ مسمتر و در مدل عددی به‌ترتیب ۰۲/۵۸ ۰۲/۲ 
٩‏ ۱/۲۶ متر محاسبه گردید که اين نتایج نشان می- 
دهد مدل عددی به‌ازای دبی‌های فوق دارای بیشترین 
رصد خطا در کمترین دبی به‌میزان ۱۱/۸۶ بوده‌است. 
نتایج مدل عددی و مدل فیزیکی باتوجه به میانگین 
1۳0/92 و ۵۷005۳ در جدول (۱) نزدیک به یکدیگر 
هستند و می‌توان به نتایج مدل عددی فلوئنت اعتماد 
کرد. در شکل‌های (4) و (1) داده‌های اندازه‌گیری مدل 
فیزیکی و محاسبات مدل عددی تغییرات عمق جریان 
در مسیر تنداب نشان داده شده‌است. برای بررسی 
هم‌خوانی نتایج مدل فیزیکی و عددی در شکل‌های (۵) 
و (۷) این داده‌ها نسبت‌به خط آرمانی (4۵0 درجه) 
ترسیم گردید و مدل برازش داده‌ها و هم‌چنین ضریب 
تعیین آن‌ها (1) محاسبه شد که نشان‌دهنده هم‌خوانی 
قابل قبول آنها می‌باشد. اين نتایج نشان که می‌دهد 
به‌ازای دبی حداقل ضریب تعیین کمتر بوده‌است. در این 
دبی. مقادیر محاسباتی با مدل عددی ارقام کمتری 
نسبت‌به مدل فیزیکی به‌دست داده‌است. هم‌چنین 
براساس معیارهای 5 و ۱5۴8 همان گونه که 
در جدول (۲) دیده می‌شود. مقادیر این ضریب برای 
کليةٌ داده‌ها مقدار قابل قبولی بوده‌است و نشان‌دهندة 
قانلنخ شاسبملال عددی مس باشد. شایان نذکر آبتته یا 
کاهش. فبی؛ میزان. خظای مدل: اقزایشین . داشته‌انست 
چراکه کالیبراسیون مدل برای دبی حداکثر صورت 
گرفته‌است. 

هم‌چنین نتایج نشان می‌دهند نوسانات عمی آب 
در کنارة دیوارة تنداب بیشتر از نوسانات در محور 
مرکزی می‌باشد. کلیدٌ اندازه‌گیری‌ها نشان می‌دهند 
دیواره‌های هدایت سرریز و تنداب از نظر طراحی 
کفایت لازم را برای عبور آب دارند. موضوع دیگر 
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مسئلةٌ جدایش جریان در محل تغییر شیب است که در مشاهده نگردید اما در آزمایش‌های با دبی‌های مختلف» 
کليةٌ آزمایش‌ها؛ نتایج نشان می‌دهند در قسمت تغییر در ابتدای شیب دوم به‌دلیل اختلاط زیاد آب و هوا 


شیب ناگهانی مدل هیچ‌گونه جدایش جریان از سطح مقداری عمق آب افزایش پیدا کرده‌است. 


جدول ۲ مقدار 1155 و ۱1845۳ نتایج عمق آب به‌دس تآمدة مدل عددی نسبت‌به اندازه‌گیری‌های مدل فیزیکی 


۷۳۳/۴۱ ۱۱۰/۸ ۱۵۱/۸ ۱۳۶/۷۳ ۳۳/۳۳ [6[6/3 
۰,۱۶۳ ۰,۶۲ ۰/۰۵۸ ۰/۰۵۹ ۰,۰۹ «متر) ۳۷6۲ 
۱/۸۶ ۱/۷۷ ۱۰/۰۸ ۱/۱۲ ۶2۹۸ ( ۱۱۳۱۷۹۲( 


8 90 72 54 36 18 0 
فاصله از تاج سرریز (متر) 


شکل ۶ تغیبرات سطح آب در مدل فیزیکی و عددی به‌ازای دبی حداکثر 
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شکل ۵ مقايسه داده‌های عمق آب در مدل فیزیکی و مدل عددی به‌ازای دبی حداکثر 


نشریه مهندسی عمران فردوسی سال سی و یکم شمارة چهار ۱۳۹۷ 


سهیل عباسی- حسین صمدی بروجنی- روحالّه فتاحی افچی- بهزاد قربابی ۱۷۳ 
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شکل ۶ تغیبرات سطح آب در مدل فیزیکی و عددی به‌ازای دبی 
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شکل ۷ مقايسة داده‌های عمق آب در مدل فیزیکی و مدل عددی 


به‌ازای دبی حداقل 


بررسی توزیع فشار 
اندازه‌گیری فشار به‌وسیلة پیزومتر انجام شد و ٩۱‏ 
موقعیت برای اندازه‌گیری پارامتر انتخاب گردید. این 
موقعیت‌ها به‌صورتی بود که ۲۱ ایستگاه در محور 
مرکزی و ۲۵ ایستگاه در محور کنارة دیواره در فاصلهٌ ۳ 
سانتی‌متری مدل فیزیکی یا ۱/۵ متری نمونه اصلی 
اندازه‌گیری شد. شکل‌های (۸) و (۱۰) فشار را برروی 
سرریز در محور مرکزی تنداب به‌ترتیب به‌ازای دبی‌های 
حداقل و حداکثر نشان می‌دهند. روند هر دو نمودار 
مشابه بوده و از تاج تا پای تنداب مقدار فشار کاهش 
پیدا کرده‌است. روند نمودارها در این شکل‌ها نشان می - 
ها که وی مه تفای ما3 ٩‏ تا ۴ بقدان فشای با 


نشریه مهندسی عمران فردوسی 


شیب تندتر کاهش پیدا کرده‌است و به‌نظر می‌رسد علت 
کاهش فشار در طول مدل فیزیکی سرریز سد سورک 
به‌صورتی بود که در هیچ‌یک از موقعیت‌های نصب 
پیزومتر» فشار منفی مشاهده نیقی باتوحه به وجود 
درز و شکاف و اهمواری به‌وجود آید. ایجاد فشار 
منفی که زمینه‌ساز کاویتاسون می‌باشد دور از انتظار 
نیست. علاوه‌بر آن هرچه دبی کاهش پیدا کرده‌است 
نوسانات مقدار فشار در طول سرریز افزايش پیدا 
کرده‌است. این نوسانات از دبی چهارم به بعد نمایان‌تر 
شده‌است. احتمال می‌رود باتوحه به این نکته که با 
کاهتن. ذبی: عمق. آب نیز کاهشی. بیدا کرده‌استه 
عملکرد پیزومترها برای اندازه‌گیری فشار در عمق کم و 
نوسانات فشار در طول سرریز در محور مرکز کمتر از 
محور کناری بوده‌است و علت این نکته می‌تواند 
علاوه‌بر ضعف پیزومتر در اندازه‌گیری فشارهای کم 
که دیواره‌های هدایت و امواج عرضی (۷۸۷ 0۲059) 
تولیدشده برروی جریان و سطح آب باشد که در 
اندازه گیری پارامتر فشار: در سرعت بالا تأثیر داشته‌اند. 
هم‌چنین برای بررسی هم خوانی نتایج مدل فیزیکی و 
عددی در شکل‌های (۱۰) و (۱۱) این داده‌ها نسبت‌به 
خط آرمانی ( 4۵ درجه) ترسیم و مدل برازش داده‌ها و 
هم‌چنین ضریب تعیین آنها (1) محاسبه شد که 
نشان‌دهندة هم‌خوانی قابل قبول آنها می‌باشد. این نتایج 
نشان می‌دهند به‌ازای دبی حداقل ضریب تعیین کمتر 
بوده‌است. در این دبی. مقادیر محاسباتی با مدل عددی 
ارقام بیشتری نسبت‌به مدل فیزیکی به‌دست داده‌است. 
هم‌چنین براساس معیارهای 5/5۲ و ۱۱۳/۲ 
همان‌گونه که در جدول (۳) دیده می‌شود مقادیر این 
ضریب برای کلیه داده‌ها مقدار قابل قبولی است و 
نشان‌دهندةه قابلیت مناسب مدل عددی می‌باشد. شایان 
داشته‌است جراکه کالیبراسیون مدل برای دبی حداکثر 
صورت گرفته‌است. 


سال سی و یکم شمارة چهار, ۱۳۹۷ 


فشار محاسبه شده 8 فشار اندازه گیری شدهحو 


1098 90 72 54 36 18 0 
فاصله از تاج سرریز (متر) 


شکل ۸ تغییرات فشار اندازه گیری‌شده و فشار محاسبه‌شده به‌ازای دبی حداکثر 
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فشار در مدل فیزیکی (متر آب) 


شکل ٩‏ مقایسه داده‌های فشار در مدل فیزیکی و عددی به‌ازای دبی حداکثر 


فشار اندازه گیری شده 72 فشار محاسبه شده 4 
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فاصله از تاج سرریز (متر) 


شکل ۱۰ تغییرات فشار اندازه گیری‌شده و فشار محاسبه‌شده به‌ازای دبی حداقل 


نشریه مهندسی عمران فردوسی سال سی و یکم شمارة چهار ۱۳۹۷ 


سهیل عباسی- حسین صمدی بروجنی- روحالّه فتاحی افچی- بهزاد قربای ۱۷۵ 
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فشار مدل فیزیکی (متر آب) 


شکل ۱۱ مقايسة داده‌های فشار در مدل فیزیکی و عددی به‌ازای دبی حداقل 


جدول ۳ مقدار 10155 و ۱1845 نتایج فشار به‌دست‌آمده در مدل عددی نسبت‌به اندازه‌گیری‌های مدل فیزیکی 


۱۵/۸ ۱۳/۸ ۷۳/۴۱ 
۰/۰۴۸ ۰/۰۵۱ ۳۵۹ 
۸/۵۵ ۱۱/۶۵ ۱۴ 


بررسی تغییرات سرعت 
ها( زر ( 0 یانش سرت روا :زا 
برروی سرریز در محور مرکزی تنداب به‌ترتیب به‌ازای 
دبی‌های حداقل و حداکثر نشان می‌دهند. روند هر دو 
نمودار مشابه بوده و از تاج تا پای تنداب مقدار سرعت 
افزایش پیدا کرده‌است. باتوجه به اين نتایج بعد از 
تغییر شیب ناگهانی تنداب» سرعت جریان با آهنگ 
تندتری افزایش افته‌است. بیشترین سرعت 
اندازه گیری‌شده جر دبی. ۲۳۹/۲ و جر فاصلة ۱۰۷۸۳ 
متری از تاج سرریز به‌دست آمده‌است و مقدار آن در 
مدل فیزیکی ۲۳/۱۰ متر بر ثانیه اندازه‌گیری شده‌است. 
این مقدار در مدل عددی فلوئنت معادل ۲۶/۵۸ متر بر 
ثانیه محاسبه گردید که نشان‌دهندهٌ خطای 7/۶ درصدی 
در مدل عددی است. وجود سرعت‌های بالا در این 
تنداب درصورتی‌که در بتن این قسمت درز و شکاف و 
ناهمواری به‌وجود آید. زمینه‌ساز وقوع کاویتاسون می- 
باشد. هم‌چنین برای بررسی هم‌خوانی نتایج مدل 
فیزیکی و عددی در شکل‌های (۱۳) و (۱۵) این داده‌ها 


سال سی و یکم» شمار چهار, ۱۳۹۷ 


۱۶/۷۳ ۳۳۴۳/۳ وان 0 
۳( ا 2 (متر) ۲۷/5۲ 
۶2۳۷ ۳۴ ۱۱۵/۱۱۰۱۸ 


برازش داده‌ها و هم‌چنین ضریب تعیین آنها (22) 
محاسبه شد که نشان‌دهندة هم‌خوانی قابل قبول آنها 
می‌باشد. این نتایج نشان می‌دهند که به‌ازای دبی حداقل 
ضریب تعیین کمتر بوده‌است. هم‌چنین براساس 
معیارهای ۲9۳ و ۱۷۷5۳ همان گونه که در جدول 
(4) دیده می‌شود. مقادیر این ضریب برای کلیهٌ داده‌ها 
مقدار قابل قبولی انستت و نشان‌دهنده قابلیت مناسب 
مدل عددی می‌باشد. شایان ذکر است با کاهش دبی. 
میزان خحطای مدل افزایش داشته چراکه کالیبراسیون مدل 

باتوجه به میانگین ضریب نش - ساتکلیف و 
میانگین درصد خطای داده‌های سرعت. در جدول (۱) 
نتایج عددی تشابه خوبی با مشاهدات آزمایشگاهی 
داشته‌است. شکل‌های (۱۱ و ۱۲) مقايسةٌ نتایج سرعت 
در مرکز و دیواره را نشان می‌دهد. به‌دلیل نوسانات 
موجود در کنارهُ دیواره مدل نتایج اندازه گیری‌شده در 


محور مرکزی دارای ضریب تعیین (187) بهتری می‌باشد. 


نشریه مهندسی عمران فردوسی 
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سرعت در مدل فیزیکی (متر بر ثاتیه) 
98 90 72 54 36 18 ۵ 0 
فاصله از تاج سرریز (متر) 


شکل ۱۵ مقايسهة داده‌های سرعت در مدل فیزیکی و عددی (دبی 
شکل ۱۲ تغییرات سرعت در مدل فیزیکی و عددی (دبی اول) 


پنجم) 
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سرعت در مدل قیزیکی (متر بر ثانیه) 


شکل ۱۳ مقایسه داده‌های سرعت در مدل فیزیک و عددی 
(دبی اول) ۱. با اندازه گیری و محاسبهةٌ فشار در ۶ دبی متفاوت 


با بررسی نتابج مدل فیزیکی و عددی می‌توان به نتایج 


زیر دست یافت: 


به‌ویژه دی طرح؛ در طول سرریز سد سورک 
هیچ گونه فشار منفی مشاهده و محاسبه نگردید. 

۲ بررسی پارامتر عمق جریان در ۲۱ موقعیت مورد نظر 
در مسیر تندابت سرریز سورک نشان می‌دهد که 
دیواره‌ها هدایت کفایت لازم را برای عبور دبی طرح 
دارند و در مکان تغییر شیب ناگهانی جدایش جریان 
وجود ندارد. 

۳ در این پژوهش در بین پارامترهای اندازه گیبری‌شده» 
نرم‌افزار فلوئنت پارامتر سرعت را با کمترین میانگین 


98 90 72 54 36 18 0۵ در صد خطا( 1/۲ در صد) و پارامتر عمق را با 
فاصله از تاج سرریز (متر) 


شکل ۱۶ تغییرات سرعت در مدل فیزیکی و عددی (دبی پنجم) شبیه‌سازی نمود. 


نشریه مهندسی عمران فردوسی سال سی و یکم شمارة چهار ۱۳۹۷ 


سهیل عباسی- حسین صمدی بروجنی- روحاله فناحی افچی- بهزاد قربالی 2 


۶ ضریب نش - ساتکلیف و ضریب تعیین محاسبه‌شده فازی ۷۵۳ به‌خوبی می‌توان پارامترهای فشار 


برای پارامترهای موردنظر در محدوده قابل قبولی 
توانایی نرم‌افزار فلوئنت را برای شبیه‌سازی جریان 
برروی سرریزهای با تغییر شیب ناگهانی اثبات می - 
کنند. باتوجه به ضریب تعیین لخادت اه برای 
داده‌های مدل عددی, مشخص گردید مدل عددی 


سرعت و عمق آب را برای بررسی تنداب‌های با 
تغییر شیب ناگهانی شبیه‌سازی نمود. هم‌چنین از 
مدل فلوئنت واسنجی‌شده در این تحقیق می‌توان 
برای بهبود عملکرد سرریز و تنداب سد سورک 
استفاده نمود. 


پارامترهای محور مرکزی تنداب را با دقت بیشتری 

سپاس گزاری 

۵ باتوجه به نتایج شبیه‌سازی‌شده و معیارهای 1۳15۳ تحقیق حاضر با حمایت معاونت پژوهشی دانشگاه 
شهرکرد به انجام رسید که بدین‌وسیله مراتب تشکر و 


قدردانی اعلام فین کو اه 


نسبت‌به داده‌های کنارهُ دیواره شبیه‌سازی نموده‌است. 


و 1180115 می‌توان ادعا نمود در نرم‌افزار فلوئنت 
بااستفاده از روش حجم محدود. مدل آشفتگی ع-عل 
0 روش تصیح فشار ٩1۳]‏ و روش دو 


مراجع 


امد ملصا فعصتمه۲ ۳۳۵۲۵07 مه ماممصمن ]۵ ممامتاتافعنا مملاهاا هن ری ۷۲۰ ری #مقاع02 ...1 
٩0.1, ۲0.24-29, )1969(,‏ و3 ,۷۵۱ رعط601صلعمظ مه زومامصطمع ]هب۳0 رقممتاتممن) 
۰ 3۳0 .1مصقن) القصصو گم صعتوعنز ررکلطونا) ممتاعصصهآهع ۵۶ بنجعتبظ تمتمنصا مطا ۵۶ اممصانهممن1 ولا ...2 
1978(۰) ,8600 ۲18۸۵ 
۳ _ مهری, م. «مدل‌سازی رفتار هیدرولیکی جریان روی سرریز سد بالارود با مقیاس کوچک» پایان‌نامة کارشناسی ارشد. دانشکدهٌ علوم و 
مهندسی آب. دانشگاه شهید چمران اهوان (۱۳۸۵). 

۶ عارف پون م.. «بررسی هیدرودینامیک جریان و منحنی تراژکتوری جت پایین‌دست در پرتابة جامی‌شکل (مطالعة موردی سد بالارود)». 
پایان‌نامة کارشناسی ارشد. دانشکدة علوم و مهندسی آب. دانشگاه شهید چمران, اهوان (۱۳۸۸). 

قر - دورف ام «زفار هیا‌رولیکی. انب برروی سرزی آوتخی شددهانه ی تست زابظه دیرت افل و بهتهسازی دیرارشای هذایت (مطالتاة 
موردی سد بالارود). پایان‌نامةٌ کارشناسی ارشد دانشکده علوم و مهندسی آب. دانشگاه شهید چمران, اهوان (۱۳۸۸). 

۰ وه [م0مصاهمتمصصار له هفرط ومسالزوگ ممو مرن ما۳۳ رصمقصومل مه ۷ رمعه927 .6 
2001(۰) ,640-649 ۲0۰ ,۱0.8 ,127 ۷۵۰ رو ۵۳۱۵۵ افیا ول 

رن 0۱۵ 0اطوتع۲۱ ماتصد گه درهسالزنگ من ما۲ تقممتاماه1 ,رل .۲ ,رلتوعمن ...7 
.(1965) ,۳(.155-173 ,0.6 ,91 ۷۵۰ بتل)۸8 ۲۱۷19109 

0 عصنام۷۵0 م۲۱ لم‌زمصتا تمممتوممصزنا ممط .1‏ رن بط رعزععاعزک1 فصه۸ .ظ نک بلط متعوان .8 
,0.775-7864( ,۱۱0.5 ,36 ,۲۷۵۱ ع0ص ۲۵6 ]0 ول ,ض9نا ,موی ردساانوگ آه ممتعتاو۳۴ 
1998(۰) 

۶ ۱۷۱۵۵6۱ لهمزمصتاا آمصمنوممصطزنا و ربا بتکلهععع] 4ص2 ,۷۲ رتحفاوظ تتصتا رآ بتطمه۳7 52 بیک رتتصقضا .9 
.(1999) ,۵369-375 ,۱۱0.4 ,125 ۷۵۱۰ ,۵۱۱6۵۲۱8۵ ع ۵ ۵۴ 0۵ ۲۱۵۱۲۰ جب 51۱ 


سال سی و یکم شمارة چهار ۱۳۹۷ نشریه مهندسی عمران فردوسی 


۱۷۸ مدل‌سازی پارامترهای جریان در مسیر تنداب سرریزهای با تغییر شیب ناگهانی... 


۱۳۵ ت۳۱ آمممتاها نموم وا دره لزوگ مین ما ۵۶ ۱۷۲۵0۵1 ,.۷ بع۲ز۵ز 20 ک .۲ رتتفهتفصصئل 
(2012) ,164-175 .۲۵ ,۱0.3 ,167 ۷۵۱ روتمعصلوص لت اه ممتاتتاعصا مطا ۵۶ فعصتلععع۳۵ظ 


10 


۱ حیدری» 1 «شبیه‌سازی جریان برروی سرریز اوجی بااستفاده از نرم‌افزار فلوئنت). پایان‌نامة کارشناسی ارشد. دانشگاه تربیت مدرس. 


تون نوی 0هممماو وحم( م۵0۵2 طا مامعله ق۰962 ۳۱ مرعمفصفط 220 .خه ر2ع[0082۵) 
۱۵۱۵۵ ۲۷۵۱۵ ۱ الا و0 0۵0۵۳۵۵ ۸۷۵۵۵09۵ ۵ 712 بر ,وسوا تعتحسستنه ما معتق‌تاو 
2004(۰) ,2112ادننظ راعومت 06010 

0۴ 0۱۲۱۵ ,5۵۱۱7۷۵۷۹۲ ۵۵۵0و ۵۲ ومتاوتمام00212 ۱۱۷۵-۰۵۳۵۵8۵ ربا ۷۷۰ ,۳12۵6۲ 200 ,۶ ,وع0ظ 
٩0.9, ۳0. 661-670, )2003(۰‏ 129 ۷۵۱ ۲۱۱۵۱۱۱6۵۳۱۸۵ 

0 0۶ ط۱2)10 نی آهع1م نار ممصمآن‌هاننع ت10 [ممجور 4تبطاع من متام ۰۲۷ ,۲۱ پتاشنا 24 .0 رون .0۵۰ بصفط) 
2002(۰) ,683-688 ۲۵۰ :128 ۷۵۱۰ رمع ]۵0 م0 0۱۵10 ره/ 50111 

٩01112۷,‏ ۳۲۱۵۷۲ و0۲ و۲۵ 2 10۲ فممتاها نی لهمزمصتار برد قمع تاد لم)ممصتهص ۲ ربظ رتطهع دنز 
2004(۰) ,899-907 .00 و132 ۷۵۱۰ رو و ۵۴۳۱۵ افیا ول 

۴ ص ومنمهطمع۱ متا آمممتاهادم‌صمی. _ ریب رتعطماعاط لعج رل باأتطعصصعه1. رین مطموتعصض 
1984(۰) ,۷۲۵۲ ۱۱۵۲ ,0۵۵۵۲20 منطو زانط وتمطام‌وتجصهم ۲ مومع 1 

۰ همم |ناط یب[ مج 0مطاع۷ متام ممتاجم2تلقصصمصف مه ما ممتامتتلنق] ریما رطف 


1993(۰) ,۰۲1۷۲-107 یبال صهمجصع]۱۷ آمعتصطاهع 1 رغصا )صمیاا۲ 
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13. 


14 


15. 


16۰ 


17۰ 


نشریهُ مهندسی عمران فردوسی سال سی و یکم شمارة چهار ۱۳۹۷ 


